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ON A CLASS OF PARTICULAR SOLUTIONS OF THE PROBLEM OF 

FOUR BODIES* 

BY 

FOREST RAY MOULTON 

Introduction. In 1772 Lagrange published a celebrated memoir on the 
Problem of Three Bodies, which contained all the solutions in which the ratios 
of the mutual distances of the bodies are constants. He found two distinct 
configurations. In the one, the three bodies always lie in a straight line ; in 
the other, they are always at the vertices of an equilateral triangle. Their 
distribution upon the line depends upon their masses, being explicitly de- 
fined by the real positive root of a certain quintic equation. The equilateral 
triangular configuration is possible for all distributions of the masses. In both 
cases the three bodies move in the same plane, in conic sections with respect to 
each other or with respect to their common center of gravity, and in such a 
manner that the law of areas is true for each body considered separately. 

No other periodic solutions of the motion of three or more bodies were dis- 
covered for more than a century, although many splendid papers appeared on 
the Problem of Three Bodies. A new impetus was given to the subject by the 
celebrated memoir of Dr. Hill on the Lunar Theory,! in which he discussed a 
new species of periodic solutions. He made the restrictions that one body 
should be infinitesimal, and that the finite bodies should describe circles around 
their center of gravity. For the purposes of the applications to the Lunar 
Theory he neglected the ratio of the distances of the moon and the sun, thus intro- 
ducing errors in the numerical results which in certain cases became important. 

Professor Daewin has more recently taken up the subject^ by a method 
not essentially different from that employed by Hill, and discussed a great many 
orbits in detail, not neglecting the solar parallax. He has considered only cases in 
which the infinitesimal body moves in the plane of the finite bodies with direct 
motion. 

To these must be added the masterful researches § of Poincae6, who has 



• Presented to the Society at the Columbus meeting, Aug. 25, 1899, under a slightly different 
title. Eeceived for publication Aujt. 25, 1899. 

tAmerican Journal of Mathematics, vol.!. 

t Acta Mathematioa, yol. 21. 

{ Lta MSthodes Nouvellea de la Mecanique Celeste, I and III. 

(Trans. Am. Math. Soc., u.) 17 
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proved the existence of several classes of infinite numbers of periodic solutions 
of the Problems of Three Bodies. His results are not given explicitly, and not 
all of the circumstances of the motion are known as they are in the previous 
cases. 

This paper treats of the Problem of Four Bodies, three of which are finite, 
moving in circles according to one or the other of the solutions of Lagrange, 
while the fourth is infinitesimal. The motion of the systems is referred to axes 
rotating with uniform angular velocity. An integral, similar to that of Jacobi 
in the Problem of Three Bodies, is found, which leads to the definition of the 
surfaces of zero relative velocity. For certain initial conditions these surfaces 
have double points, all of which are found. 

It is shown that if the infinitesimal body be started at any double point of the 
surfaces, with zero relative velocity, it will always remain in the same position 
with respect to the finite bodies. Therefore they are particular solutions of the 
motion of the system, and it is true that they are the only ones in which the mu- 
tual distances of the bodies are constants. There are eighteen solutions of arbi- 
trary periods of revolution in which the finite bodies lie in a line, and ten in which 
they are at the vertices of an equilateral triangle. 

§ 1. The differential equations of motion. 

Suppose the masses of the finite bodies are m^, m^, and m^ respectively. Sup- 
pose the origin is taken at the center of gravity of the system. Let the coordi- 
nates of m^ be ^^, ^^^, ^., i = 1, 2, 3, and of the infinitesimal body, 6, tj, and ^. 
Suppose the units are chosen so that the Gaussian constant is equal to unity. 
Let 



and 

f'l P2 Pi 

Then the differential equations of motion of the infinitesimal body are 

f d'^ _ dU 



(1) 



d^i _dv 

de ~ d-q ' 

dX _ dU 
df ~ ~"d^ ' 

The motion of the system is referred to ccy-axes rotating in the f ij-plane by the 
substitution 
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(2) 



^, = X. COS d — y. sin d, 
yj. = X. sin d + y^ cos d, 



where ^ is a function of the time. 

Equation (1) is transformed by this substitution into 



W 



(8) 



dd_d^ 
dt dt 



(dOV 



d^ 



dU 

y=-dx^ 



W ^ dt dt~ \dt) y'^ df"^- dy' 



dh 
d? 



dz ' 



If the axes rotate with the uniform angular velocity, n , the system (3) 
becomes 



(4) 



d'x 



dy 



dU 



dv dt dx 



m^{x — x^ 



f\ 



mix — x^ mix — X3) 



P\ 



de ^^""dt -"'y+ dy-'^y p\ 



'fnly - 2/2) '"hiy - 2/3) 



p\ 



p\ 



dh 
df 

The substitution 



d U mlz — Sj) mlz — %^ mlz — a,) 



(>\ 



Pi 



(5) 



dx , 

57-^*2' = ^' 

dy 
dt 

dz 
dt 



+ nx = y', 



H= \ (x' + ny)' + iiy' — nxf + z'^^— in\x^-\- y^) — U 
transforms the system (4) into the canonical system 
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(6) 
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dx dH 


dx' dH 




dt ~ dx' ' 


dt dx ' 




dy dH 

dt ~ dy' ' 


dy' dH 

dt ~ dy ' 




dz dH 

dt ~ dz' ' 


dz' dH 

dt ~ dz ' 



§ 2. Integrals of the differential equations. 

Suppose the angular velocity of revolution of the finite bodies, which we now 
assume move in circles according to one or the other of the Lagrangian solutions, 
is also n ; then x^, y^, z. are constants. Therefore, from system (6) there results 
the integral ^ + C = 0, or, in full, 



/-TN (^^y (^y\ (^^y 2/ 2 2x 2mj 2m. 



^'_a 



Let r^ = y/ix — x^^ + (y — Vif ! then, if we suppose that the infinitesimal body 
moves in the a;2/-plane, the system of equations ( 6 ) reduces to the fourth order, 
and the integral becomes 



.„, / dxV /' dyY _„ ,^ 2m. 2m. 2m. 



- a 



If one integral more could be found under this last hypothesis, JaCOBI's theory 
of the last multijMer would be applicable, and the problem could be completely 
solved. However, this system of differential equations is of the type considered 
by Bruns in Acta Mathematica, vol. 11, for which he proved that there are 
no new algebraic integrals. It is equally of the type considered by PoiNCARft 
in Les MetJiodes Nouvelles de la Mecanique Celeste, vol. 1, chap. V, for 
which he proved that there are no new uniform transcendental integrals, even 
for very small relative masses of two of the three finite bodies. 

§ 3. The Lagrangian solutions when the bodies revolve in circles. 

Suppose first that the three bodies move so as always to lie in a straight line. 
Suppose the order of the bodies is wi, , m^ , m^ . Suppose the distances between 
m^ and m^, m^ and m^, m^ and m^ are respectively 1, a, and a^. Therefore 
Oj = 1 + ttj . Then the equation for the determination of a^ is * 

(9) (tMj + m^)a\ + (3m, + 2m^a\ + (8m, + m.^al — (m.^ + Smj)^, 

— (2m2 + 3m3)a, — (m^ + mj = 0., 



* Laplace, Mec. CSl., vol. V, p. 310, or Tisseeand, Mec. CM., vol. I, p. 155. 
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The a^ to be chosen is the single real positive root of (9). Then the coordi- 
nates of the three bodies are given without ambiguity by 



(10) 



, cc. ^= x^ 



*^2 

JU,, ^^ JCn —~ (ai 



Vi 



1» 

yi = yi = 0- 



As there are three arbitrary units, of which only two have been thus far se- 
lected, the angular velocity may be chosen arbitrarily, and the unit of mass then 
determined. Or, if the unit of mass is taken arbitrarily, the angular velocity 
is given by the equation 



(11) 



?r = mj + wij -f m^ 



\ «2 <A y 



There are two other different arrangements of the three bodies upon a straight 
line, in which their order is cyclically permuted. Their coordinates are de- 
termined in a manner similar to that employed above. 

In the second Lagrangian solution the three bodies always lie at the vertices 
of an equilateral triangle. A side of the triangle being taken as the unit of 
length, the expression for the angular velocity is given by 



(12) 



m^-\- m^-\- my 



(13) 



Suppose that the axes are so chosen at the origin of time that m^ is on the 
a>axis ; then the coordinates of the three bodies are determined by the equations 

2/i = 0> 

'^{Ux + ™22/2 + '^Hi = 0' 

(a,^ - x^^ -f (y, - yy^ = 1, 
(a,3 _ x^^ + (^3 - y,)2 = 1, 
(a,^ _ a;,)^ + (y^ _ 2/3)2 = 1. 

There is evidently but one distinct configuration of this type for three given 
bodies. 

§4. Surfaces of zero relative velocity in the case of the straight-line Lagran- 
gian solutions. 

The relative velocity of the infinitesimal body being represented by V, equa- 
tion (7) may be written in this case 
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(14) V'^nXx'' + f) + ^"''1 -L ^"^-^ 



\{x-xj + y'+zr^' • l{x-x,y + y'+z'r' 

, ^ c 

This equation serves to give the relative velocity of the infinitesimal body at 
every point of rotating space ; or, if the velocity be chosen arbitrarily, it gives 
the surfaces at which the infinitesimal body will move with the given velocity. 
The most interesting surfaces are those at which the velocity is zero, for they 
separate those portions of relative space in which the infinitesimal body may 
move from those in which it cannot. Setting V equal to zero we have 

(14) nXx" + f) + j(^_^j2^'^2^,^2|K + i^-z.-^r:^rf:f^\^ 

-^ {(x-x,y + y'-\-z'\^^^- 

The surfaces of which (14)' is the equation are symmetrical with respect to 
the xy- and cca-planes, and a line parallel to the g-axis will pierce them in two (or 
no) points. Moreover, the surfaces are contained within a cylinder whose axis 
is the a-axis, and whose radius is \^ C/n, to which certain of the folds are 
asymptotic at a^ = oo. 

The intersection of the surfaces with the xy-plane is given by 



2toj 2ni^ 2 



m. 



(15) nV + f) + "P + =P + ":r-^ = C'. 



'2 '3 



For large values of x and y which satisfy this equation we have approximately 

n\x' + y") = G. 

Therefore for sufficiently large values of C one branch of the curves in the 
ajy-plane is approximately a circle near the asymptotic cylinder. 
For small values of x and y satisfying (15) we have 

m, m„ »i„ C 

L _L — £- _[_ — 2- =_ 

A. A. A. V 

'l '2 '3 

This is the equation of the well-known equipotential curves for three centers 
of force, and for sufficiently large values of C they consist of closed ovals 
around each of the bodies. For smaller values of C these ovals coalesce by 
enlarging, while the large oval mentioned above diminishes in size. 

The approximate equations of the curves in the X2-plane, for large and small 
values respectively of x and z are 

mV = c, 

wij m^ w, C 

|(a,_ccj^ + »2p + |(jg_a,^2^g2p + ]{x-x^Y + z^^' ^ "2 " 
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The first is the equation of two straight lines parallel to the g-axis ; the sec- 
ond is the equation of equipotential curves. Therefore the locus in the cca-plane 
for sufficiently large values of C is composed of curves asymptotic to the inter- 
section of the asymptotic cylinder with the xa-plane and of ovals around each of 
the three bodies. 

The approximate equations of the curves in j/a-plane for large and small 
values respectively of the coordinates are 

nY = a 



"T /-„2 1 ..2 1 Ji,\M r 



(x^ -f 2/^ -f z'Y' ^ {x\ + y' + z'f -^ {xl + f + zY' ~ 2 

The first is the equation of two straight lines parallel to the g-axis. If the 
bodies lie in the order m^ , m^ , m^, then x^ will be much smaller than x, and 
Xj , unless there is a very great difference in the masses m^ and m^. If wt 
and TWj are nearly equal the last equation becomes approximately 

m.^ C 

which is the equation of a circle. If m^ is very much greater than m^ , then 
aj, is small compared to x^ and x^ , and we have approximately 

m, C 



{x] + y' + z'f ~ 2 

which is also the equation of a circle. 

Thus we conclude that for sufficiently large values of C, the surfaces are com- 
posed of a closed fold, roughly spherical in form, around each of the finite 
bodies, and curtains hanging from the asymptotic cylinder symmetrically with 
respect to the .ry-plane ; for smaller values of C, these folds and curtains 
unite ; and for sufficiently small values of C the surfaces are composed of two 
parts symmetrical with respect to the xy-plane but not intersecting it. 

The motion of the infinitesimal body is real within the small folds around the 
finite bodies, and outside of the curtain hanging from the asymptotic cylinder, 
while it is imaginary in all the other portions of relative space. Therefore, if 
the infinitesimal body is within one of these closed regions at the origin of time, 
it will forever remain there. 

Points on the surfaces may actually be computed most readily by first deter- 
mining the curves in the xy-plane, and then finding, by methods of approxima- 
tion, values of z' which will satisfy (14)'. This is true since approximate values 
of z^ are known from the general nature of the surfaces, when the curves in the 
ajy-plane are given. 
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If (15) is rationalized and cleared of fractions the resulting equation is of the 
40th degree in the variables x and y, and may be used for the determination of 
points on the curves. The difficulty of locating the roots of an equation of so 
high degree is supplemented by the presence of foreign solutions which have been 
introduced in the process of rationalization. The latter difficulty can be avoided 
entirely and the degree of the equation notably reduced by transforming to bi- 
polar coordinates. 

If Wj is supposed to lie in the negative direction from the origin on the xaxis, 
it is found without difficulty that 

(16) a-- + 2/^ = ^^i — x^rl + x^x, . 
Then (15) becomes 

(17) n^ix,i - x^ + 1;'^' + ^ + -^ = C - r^x^, . 

'i '2 '3 

r^ can be expressed rationally in terms of r^ and r^ by an equation of the second 
degree. Therefore equation (17) may be written as a polynomial in r^ and r^ of 
the 10th degree. Points on the curves may be found by assuming convenient 
values of one of the variables and finding the real positive roots of the resulting 
equation, which will be of the 8 th degree in the other variable. The points of 
intersection of the circles around niy and w^ as centers with the assumed and 
computed values respectively of r^ and r^ as radii, will be points of the locus. 

§ 5. Double points of the surfaces. 

It follows from the symmetrical forms of the surfaces that, as they degener- 
ate and coalesce with decreasing values of (7, all of the double points will ap- 
pear in the aj^z-plane. Therefore the double points of the curves in the a;y-plane 
alone may be sought. That is, the values of C for which double points occur 
and their positions must be determined from the equation 

ri8^ nHx'' + v')+ '^^ + --' 

+ l{x-x^f+y^^'-^- 



The conditions for a double point are 

dC _ 2 m,^(x — x^) m.j^x — x^ mj^x — x. 



(19) 



dx '" r'l r\ r* 

dy^^y r\ r; rl 



= 0. 
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There are certain double points on the .r-axis. They are determined by the 
equation 

(20) n?x - "^i( ^ - ■'^i) _ "^l ^ - ^ 2) ' "h{^ - ^3) ^ 0. 



K*-a^i)'5^ {(3^-2^2)'^ K^'' 



a;. 



The left side of this expression, considered as a function of x, changes sign 
once, and only once, in each of the intervals + <» > a: > .T3 , ajj > a; > ajj , aj^ > a; 
> «j , a!j > a; > — 00 . Therefore there are four double points of the curves 
on the a;-axis, one in each of the intervals mentioned. Expressing equation (20) 
in terms of r^ as the variable and clearing of fractions, we have a rational equa- 
tion of the 7th degree in r^ It has but one real positive root for the form which 
it takes in each of the intervals ; therefore it locates the positions of the double 
points, which are on the a;-axis, uniquely. By substituting the coordinates of 
these various points in (18) the corresponding values of C are determined. 

The double points for which y is different from zero are found by solving (19) 
after dividing the second equation by y. By simple combinations of the two 
we find 



(21) 



U 3 ^^ _ ^^)2 + y2|3 j-^^ _ -^.^^2 + 2/^ Ji - ' 

^'^ "" ^'^ + l{x - x^y + y^i + {(x- x^y + fli = ' 
^^ ■' l(x-x,y + y^i+l{x-x^y + f[i 



Any one of these three equations may be derived from the other two, and the 
simultaneous solution of any two of them will furnish all of the double points 
not on the .T-axis. 

The equations (19) are equally the necessary conditions for a minimum value 
of C. From (18) it is evident that C has a minimum value, and it is easy to 
show that the sufficient conditions have not been fulfilled in the case of any of 
the solutions for which y is zero, hence it is inferred that there is at least one 
solution of the system (21). 

Suppose mj > m^ ; then x^ < 0, a;^ > ^ 3,nd a;3> a;^ > 0. The equation (c) is 
not fulfilled for any value of y^ when x is greater than x^ ; hence there is no so- 
lution in the interval x^-^ x <C + 00 . Likewise (b) is not fulfilled for any value 
of y^ when x is less than .^•, ; hence there is no solution in the interval 
a?! s X > — 00 . Therefore the solutions must all lie in the interval 
a^i < « < ^2- ^^^ small values of y', (h) and (c) are fulfilled by values of 
X near a;^ and x.^ respectively. For increasing values of y^ they are fulfilled 
by values of x continuously increasing in distance from a;, and x^ respectively; 
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therefore there is one, and only one solution in this interval. The sufficient 
conditions for a minimum value of C are fulfilled in this case, and the curves in 
the cc!/-plane have degenerated to two points symmetrically situated with respect 
to the cc-axis. There are, therefore, for one distribution of the finite bodies on 
a line, four double points on the a>axis and two situated symmetrically with 
respect to it. 

In the special case in which m^ = m^ the solution occurs for x == cc^ = 0, the 
y coordinate being given by the equations 



(22) 



2 "^l , 



The y coordinate has a maximum value for m^ = 0, and decreases as m^ in- 
creases. 

§6. Particular solutions of the motion of the infinitesimal hody when the 
finite bodies move in circles according to the straight-line Lagrangian solu- 
tion. 

The surfaces of zero relative velocity discussed in §4 are surfaces which the 
infinitesimal body cannot cross. Besides, it is easy to show that the lines of 
effective force are orthogonal to these surfaces. It follows that if the infinites- 
imal body is placed at a point in a surface of zero velocity it will move at first 
orthogonally to the surface. At the double points it seems that there would 
not be a tendency to move in any direction and it might be inferred that the 
body would remain relatively at rest. 

The direct proof of this is that if the coordinates of the infinitesimal body 
fulfil the conditions for double points, (19), they also identically verify the 
differential equations of motion, (4). 

There are thus, for each of the three permutations of the order of the finite 
bodies, four solutions of the problem in which the four bodies all lie in a straight 
ine, and two in which only the three finite bodies lie in a straight line. Therefore, 
for three given finite bodies and one of infinitesimal mass, there are eighteen, 
and no more, different solutions of arbitrary period of revolution, such that the 
finite bodies lie in a line while the mutual distances of the four bodies are con- 
stants. 

§ 7. Surfaces of zet'o relative velocity in the equilateral triangular La- 
grangian solution. 

Suppose the three finite bodies revolve in circles so as to remain at the 
vertices of an equilateral triangle, and that the axes are chosen so that the x-axis 
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passes through m^. Then the equation of the surfaces of zero relative velocity 
is found from (7) to be 



2/ 2 2\ 2m, 2m 



(23) 



l{x - x,f + f + zY' ^ \{x - x,f + {y- y,f + z^' 



The surfaces are symmetrical with respect to the aiy-plane, and are enclosed 
within a cylinder whose axis is the »-axis and whose radius is n/ C/n, and to which 
certain folds are asymptotic at »^ = oo. 

The equation of the trace of the surfaces upon the ajy-plane is 

J , J 2. 2m, 2m2 

'"^"^ +y)+l{x- x,f + y^ + {{x - x,y + {y- y,f]y' 

+ ^(x-x,f+{y-y,f\^ ~ 

For respectively large and small values of x and y satisfying this equation 
we have approximately 

n'ix' + f) = C, 
2m, 2m2 2m3 



\{x - x,y + yT' ^ \{x- x,f + {y- y^^' ^ \{x - x,f + {y - y,f 

The first equation shows that, for sufficiently large values of C, one branch of 
the curves is approximately a circle lying near the asymptotic cylinder. The 
second is the equation of the equipotential curves, which, for sufficiently large 
values of C, consist of closed ovals around each of the three finite bodies. For 
smaller values of C these smaller ovals coalesce by enlarging, while the large oval 
unites with the smaller ones by decreasing in size. 

The trace of the surfaces in the CK^-plane is composed of two curves asymptotic 
to the lines of intersection of the plane and the asymptotic cylinder, together with 
an oval around m,. If the a;- axis had been directed so as to pass through either 
of the other two finite bodies similar results would have been true. Therefore, 
we conclude that, for sufficiently large values of C, the surfaces are composed 
of curtains hanging from the asymptotic cylinder, and closed folds approximately 
spherical in form around each of the finite bodies. 

The motion of the infinitesimal body is real within the closed folds and out- 
side of the curtains. If the infinitesimal body were within one of these closed 
regions at the origin of time it would forever remain there. 

In order actually to compute points on the surfaces, it is advantageous to 
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transform the equation of the curves in the a3(/-plane to bipolar coordinates. 
Points on the surfaces can then be rapidly found by methods of approximation. 

§ 8. Double points of the surfaces. 

It follows from the symmetry of the surfaces that the double points which 
appear as C decreases will all lie in the a;?/-plane. They will evidently appear, 
(a) where the small ovals around the finite bodies first touch each other ; (6) 
where the small ovals first touch the large oval ; (c) where the small triangular 
curves, whose sides are formed from two of the small ovals and the large one, 
vanish from the xy-plane ; and (d) where the triangular curves, whose sides are 
formed from the three small ovals when they unite, vanish from the icy-plane. 
The points (c) and {d) occur for minima of C as defined by (24), the absolute 
minimum depending upon the distribution of the masses in the finite bodies. 

The conditions for double points are found from (2-4) to be 



(25) 



OX r, rl rl 



^^-„,2,. '^iV '"\^y-y^ ^ly-Vi) 



= n^y — 



dy r'j r^ 



= 0. 



In order to exhibit most clearly the distribution of the double points and omit- 
ting to give the direct proofs from the equations, which would be similar to those 
given in the former case, let us solve numerically the simple example in which 
TO^ = TWj = Wj = 1. Then «^ = 3. Since the figure is symmetrical with respect 
to the three axes passing through the center of gravity and the three bodies 
respectively, and since all the double points lie on these axes, it will be suf- 
ficient to locate the double points on one of them. In this case the second of 
(25) is fulfilled for y = 0, and the first becomes, a side of the triangle being 
taken as unit of length. 



K''-;>r)r K"+23f)' + 4 



The real solutions of this equation are approximately x = 1.22, a; = 0, a; = — .24 
and a; = — .93. The first belongs to the class (6) defined at the beginning of this 
section, the second to the class (J), the third to the class (a), and the fourth 
to the class (c). There are corresponding solutions on the other two axes, and 
these are the only real solutions of the equations. As the point a; = has been 
counted three times it is seen that there are ten double points of the surfaces 
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of zero relative velocity in the case of the triangular solution of the problem of 
three bodies. 

§ 9. Particular solutions of the motion of the infinitesimal hody in case 
the finite bodies move in circles according to the equilateral triangular La- 
grangian solution. 

By general considerations similar to those advanced in § 6 we conclude 
that the double points are points of zero effective force, and that if the in- 
finitesimal body is placed at one of them with zero relative velocity it will always 
remain in the same relative position. 

The direct proof that these are particular solutions of the problem is that if 
the relative coordinates of the infinitesimal body are supposed to be constants 
fulfilling the conditions for double points (25), they also verify identically the 
differential equations of motion (4). 

Therefore, three finite bodies of arbitrary masses and one infinitesimal 
body may be started in twenty-eight, and no more, difi'erent ways so that they 
will revolve in an arbitrary period and preserve their mutual distances. 

Thb University of Chicago. 



